Helicoidal Fields and Spin Polarized Currents in CNT-DNA Hybrids 



G. S. Diniz/'B A. Latge,2 and S. E. Ulloa^ 

^Department of Physics and Astronomy and Nanoscale and Quantum 
Phenomena Institute, Ohio University, Athens, OH 45701, USA 
^ Instituto de Fisica, Universidade Federal Fluminense, Niteroi-RJ 24210-340, Brazil 

(Dated: December 20, 2011) 

We report on theoretical studies of electronic transport in the archetypical molecular hybrid 
formed by DNA wrapped around single- walled carbon nanotubes (CNTs). Using a Green's function 
formalism in a vr-orbital tight-binding representation, we investigate the role that spin-orbit inter- 
actions play on the GNT in the case of the helicoidal electric field induced by the polar nature of 
the adsorbed DNA molecule. We find that spin polarization of the current can take place in the 
absence of magnetic fields, depending strongly on the direction of the wrapping and length of the 
helicoidal field. These findings open new routes for using GNTs in spintronic devices. 



Carbon nanotubes continue to attract a great deal of 
attention due to their potential application in a wide vari- 
ety of electronic devices, including molecule detectors [l| 
and spin-based electronics [H, to name but a few. Ad- 
vances in synthesis allow production of GNTs with spe- 
cific metallic or semiconductor behavior, depending on 
the tube chirality, and further tailoring of their electronic 
and optical characteristics via chermcal functionalization 
with specific atoms and molecules [3]. 

Many efforts in a variety of carbon nanostructures have 
been directed to study the effects of spin-orbit inter- 
action 0-10]. Although this effect is expected to be 
weak, due to the low atomic number of carbon, it has 
been demonstrated in a variety of systems that depend- 
ing on substrates or adsorbed atoms, the spin-orbit in- 
teraction (SOI) can be effectively enhanced in carbon 
nanostructures [ll|, |T2[. In GNTs in particular, the cur- 
vature arising from their cylindrical geometry can en- 
hance SOI via orbital hybridization among neighboring 
atoms jsl, m, [lo|, 13 1 . The application of external elec- 
tric fields gives rise to Rashba SOI (RSOI) effects that 
may even create electronic 'helical' states for fields ap- 
plied transversely to the GNT [lo|. Most significantly, 
recent experiments have demonstrated that SOI can be 
quite important in determining the electronic structure 
and transport properties of GNTs [l^ . 

Ghemical functionalization of GNTs is an important 
route to manipulate their properties 15 -l^- Prominent 
among these is the hybrid formed by DNA molecules ad- 
sorbed on the exterior walls of the GNT, and found to 
arrange in a wrapped-around fashion in clock- and coun- 
terclockwise helices with pitch (coiling period) which cor- 
relates with GNT chirality [l|, [l8|-[20| . The DNA molecule 
is stabilized on the GNT by van der Waals forces, driven 
by the optimal overlap of 7r-orbitals of the DNA bases 
and of the carbon hexagons on the cylindrical nanotube 
(i.e., optimal tt-tt stacking), as well as by electrostatic 
interactions of the phosphate backbone [l|, [l9|, |20| . In- 
terestingly, as the backbone is charged, the wrapping of 
the molecule generates sizable helicoidal electric fields on 
the GNT walls. Helically modulated potentials have been 



shown to affect the electronic spectrum of GNTs 21, 22[, 
however, the role of the strong RSOI effects generated by 
these fields has not been considered. 

SOI and helicoidal fields have been identified in re- 
cent experiments to produce spin selective transport of 
photoelectrons on DNA molecules on a gold substrate 
[23|. Theoretical studies explain the high polarization 
as arising from the RSOI induced by the helicoidal field 
built in the chiral molecule, as well as by the inherent 
low- mobility of carriers in DNA [2^, in agreement with 
direct transport measurements on similar systems (iBj . 

In this letter, we present results of electronic transport 
in GNTs in the presence of RSOI originating from the 
electric field induced by the DNA molecules. We find 
that as the typical chiral wrapping of the molecules pro- 
duces a helicoidal field, this results in spin polarization of 
the current that increases quadratically with the number 
of wrappings for experimental estimates of the RSOI field 
strength. This polarization in the absence of magnetic 
fields suggest the possible applications of GNT-DNA hy- 
brids as spin polarization devices, as well as the plausible 
detection of adsorbed molecules on GNT by means of po- 
larized electrical conductance measurements. 

The GNT-DNA system is modeled using a 7r-orbital 
orthogonal tight-binding Hamiltonian that includes the 
effect of RSOI due to the helical DNA molecule [26|, 
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where c]^ creates an electron at site i with spin a (=t, i) 
and Cia- is the on-site energy, affected by the strong local 
fields. Oaa' = t^crcr' + ^^^^ {^ij • s)crcr' , iucludcs the near- 
est neighbor tt-tt hopping, t ~ — 2.9eV, and the Rashba 
coupling that depends on the electric field strength 
at the nearest neighbor atom pair 5 are the Pauli 

matrices, and Uij = z x ffij^ where z is the unit vector 
perpendicular to the GNT's surface and ffij is the vector 
connecting (z, j) [26|, [13 • For GNTs with small radius 
there is an additional SOI term arising from the curva- 
ture [gI, [io| ; in this work we focus solely on the stronger 
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RSOI produced by the electric field induced by the DNA 
molecules. 

We model as a constant Vr for those CNT atoms 
in close proximity to (underneath) the DNA helix, ef- 
fectively cutting off the long-range of the helicoidal field 
distribution. This short-range approximation mimics the 
screening of the field by the CNT and background sub- 
strate, and it does not qualitatively affect our main re- 
sults and conclusions. The electrostatic potential also 
lifts the on-site degeneracies via e^cr in Eq. (dl [io|,[28|. 




FIG. 1. (Color online) CNT-DNA hybrid: DNA wraps around 
the CNT in a helical fashion. The bottom figure shows an "un- 
zipped" version of the CNT system; the central region with 
adsorbed DNA is connected to left/right semi-infinite CNTs. 
DNA helix adsorbs with pitch le, for a given wrapping angle 
with respect to the tube axis; wrapping can be clockwise (+) 
or counterclockwise (— ). = 0,±7r/6, and ±7r/3 are shown. 

The spin-resolved conductance is calculated using a 
surface Green's function approach in real space The 
CNT device is divided into three regions: left lead, cen- 
tral conductor and right lead, as shown schematically in 
Fig. [H The central conductor contains the CNT-DNA 
hybrid, the only region under the influence of RSOI ef- 
fects; it is also connected to semi-infinite leads by nearest- 
neighbor hopping with perfect atomic matching, to sim- 
ulate the DNA adsorbed on a perfect CNT. The Green's 
function of the central conductor (omitting the spin in- 
dices) is then 

g^^'-iE) = {io± -Hc-^L- Sfl)"' , (2) 

where a/r denotes the advanced/retarded Green's func- 
tion (with energy uj± = E ± if], respectively, ^ 0), 
and E is the energy of the injected electron (the Fermi 
energy at a given doping). He stands for the Hamilto- 
nian in the central region and S/^/^ are the self-energies 

for the right/left leads, T^i = Hj^giHic, where gi is the 
Green's function for the I = R semi-infinite lead, ob- 
tained through an iterative procedure [29], and Hic cou- 
ple each lead to the central region. The spin resolved 
conductance through the entire system is given by, 

= GoTr [r^,gh,,^,T^,g^c,a'.] , (3) 



where the trace runs through the lattice sites. Go = e'^/h 
is the quantum of conductance per spin, and FJ^ are the 
couplings for the leads, related to the spin-diagonal self- 
energies by F^ = z m - Sf] [igj. 

We focus on the results for the electronic conductance 
of a metallic zigzag (9,0) CNT. These illustrate the gen- 
eral behavior for metallic zigzag tubes, sharing the linear 
dispersion close to the neutrality point (the Fermi level 
of undoped CNTs), and similar spin conductances close 
to the Fermi energy. For higher (or lower) Fermi en- 
ergy values (n- or p-doping), the secondary energy shells 
naturally contribute to the spin conductances for both 
metallic and semiconductor CNTs, although this requires 
substantial doping; we focus here on the low-doping be- 
havior. 

The wrapping angle 6 of the helical molecule around 
the CNT is defined so that 6 = corresponds to a 
molecule adsorbed parallel to the CNT long axis. In con- 
trast, a helix with = ±7r/6 or ±7r/3 follows the max- 
ima of the orbital wave function on the CNT's surface 
and produces the expected optimal 7r-stacking As 
shown below, 6 strongly affects the electronic transport 
properties on a given nanotube. 
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FIG. 2. (Color online) G^a' for CNT-DNA hybrid with zigzag 
(9,0) CNT for adsorbed DNA parallel to the CNT axis {0 = 
0); length of central region is 5.12nm. (a) Rashba parameter 
Vr = O.lt; (b) Vr = 0.2t. Solid lines are for = 0; dashed 
lines Cia = O.lt. Line colors specify G^a' in the legend; in all 
cases G^^ = G^ and G^i = G^^ for this angle ^ = 0. 

Fig. [2] shows the conductance spin components for a 
CNT with adsorbed DNA parallel to the tube axis, i.e. 
^ = 0. The length of the central tube is 4/51^7^/3 = 
5.12nm. Each panel shows results for different Rashba 
parameter Vr^ with (solid) and without (dashed) on-site 
energy modulation. One can identify important features: 
(i) The diagonal (non-spin-flip) conductance components 
obey G^^ = G^^, while off-diagonal (spin- flip) compo- 
nents are G^^ = G^^; these equalities are expected from 
general considerations for time-reversal and parity pre- 
serving fields, (ii) The on-site energy modulation breaks 
electron-hole symmetry, as the induced electric field acts 
on the affected sites as donor /acceptor local fields which 
scatter charge carriers, (iii) As Vr increases, there is a 
drastic enhancement of the spin- flip channels, as the non- 
spin- flip conductance components decrease, as one could 
intuitively expect. However, notice that symmetry (i) in 
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this case guarantees there is no spin-polarization of the 
current, regardless of the value of Vr or the length of 
DNA strand. 
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FIG. 3. (Color online) G^a' for (9,0) CNT and DNA wrapped 
4 times with Q — zLtt/G, tia — O.lt. (a)-(b) Results for Vr — 
O.lt. Notice G^^ exchanges with G^^ and G^^ with G^^ when 
wrapping direction reverses from (a) to (b). In (c) Vr = 0.2t, 
6 — +7r/6. Colors follow legend in (a), (d)-(e) Gcra' for 
tia = and = +7r/6 as function of Vr for E = 0.05t and 
E — 0.2t, respectively. 



Let us now consider the helicoidal field produced by a 
helical wrapping of the molecule. Fig. [3] shows the con- 
ductance spin components when the DNA is wrapped 4 
times with 9 = ±7r/6; each DNA loop extends le=7r/6 — 
3.84nm along the tube direction. Panels (a) and (b) show 
results for the two directions of wrapping, with Vr = O.lt; 
(c) shows results for Vr = 0.2t for the direction, (d) 
and (e) show the oscillatory behavior of the spin com- 
ponents as function of Vr for two values of the Fermi 
energy, E = 0.05t and 0.2t. A few important points are 
illustrated by these graphs: (i) There is a clear difference 
between the spin resolved components for both non-spin- 
fiip and spin- flip conductances, such that G^^ ^ G^^, for 
example. This difference is strongly energy dependent 
and can be significant for all components. The splitting 
of the curves is clearly caused by the chiral parity sym- 
metry breaking of the helicoidal field generated by the 
wrapped DNA molecule, (ii) The asymmetry exhibited 
by Gcra' is reversed when the direction of wrapping re- 
verses. We stress that this asymmetry in G^a' does not 
violate time-reversal invar iance, as the inversion under 
wrapping reversal illustrates. In other words, reversing 
the direction of the current for a given wrapping results 
in the reversal of the various G spin components, as ex- 
pected from symmetry. An important result of this heli- 
coidal field and RSOI is the possibility of spin polariza- 
tion of the current for a given bias direction, as we will 
see below, (iii) The asymmetry in G^a' is not the same 



for electrons and holes due to the on-site energy mod- 
ulation considered-this can be easily verified by setting 
^ia =0- (iv) The different spin components of G depend 
on an oscillatory manner on the value of Vr, as shown in 
panels (d) and (e), for a given length and wrapping angle 
of the helicoidal field. This behavior is reminiscent of the 
spin field effect transistor and has a similar source (soj . 
as the spin precesses as it propagates in the presence of 
the Rashba field, acquiring a net phase that is propor- 
tional to VrL, where L is the number of DNA wrappings. 
Other wrapping angles show similar behavior. 

Results on spin asymmetry and general qualitative de- 
pendence on Vr and wrapping length are similar in other 
metallic tubes (such as (12,0)), with slightly different en- 
ergy dependence. The effects on semiconducting tubes, 
such as (10,0) only appear, of course, at higher energies 
(beyond the gap). Similarly, the values and symmetries 
among the different spin components of the conductance 
depend strongly on the wrapping angle. 

We now compute the normalized conductance polar- 
ization P = ^^{G^a - GicT)/Y.^^>GcTCT'' Results for 
= ±7r/6 are shown in Fig. |4] as function of energy. 
The polarization is reversed as the wrapping direction 
is inverted, i.e. as it changes from (solid lines) to 

" (dashed)-a consequence of time-reversal symmetry 
in the problem. Notice also that the perfect asymmetry 
reversal for electrons and holes P{E) = —P{—E) is bro- 
ken when an on-site energy modulation is present, see 
inset. The polarization amplitude is strongly dependent 
on the number of wrappings, growing by nearly one or- 
der of magnitude from one full loop (black lines) to four 
full wrappings (blue). P also increases with Vr, nearly 
quadratically in this range, as shown from panel (a) to 
(b) (notice different vertical scales). 




FIG. 4. (Color online) Normalized conductance polarization 
P of a (9,0) CNT as function of injected electron energy for 
different DNA lengths with = ±7r/6. (a) Vr = O.lt; (b) Vr 
= 0.2t. Solid (dashed) lines represent results for "+" ("— ") 
direction; P increases with number of wrappings (1 through 
4 shown with different color lines). Central inset shows P for 
eia = 0, Vr = 0.2t, and = zLtt/G, displaying full electron- 
hole asymmetry, P{E) — —P(—E). 

Fig. [5] shows P as the number of wrappings L and 
Vr change, for an energy close to the Fermi energy, 
E = 0.03t. DNA parallel to the CNT's axis (dot-dashed 
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line) produces no polarization, regardless of the RSOI 
strength and length; i.e., a non-helicoidal field produces 
no spin polarization. On the other hand, non-zero polar- 
ization is evident for helicoidal fields. The polarization is 
found to oscillate for each configuration, as shown in the 
inset of Fig. [5l the peak polarization value is c:^ 35% for 
L = 8 for the system with = — 7r/3 and Vr = 0.2t. Sim- 
ilar non-monotonic behavior is found in P for a given 
and DNA length, as Vr is varied, reaching a peak value 
similar to that in Fig. [5] (not shown). This oscillatory 
behavior is in agreement with the expected phase accu- 
mulated by the spin precession, and proportional to VrL^ 
as described above. We find that the polarization follows 
an approximate relation, P :^ Pm sin^(yHL/t), with Pm 
depending on (typically Pm ^ 30% to 40%). This re- 
lation allows one to estimate the polarization for a given 
CNT-DNA hybrid, once Vr is known. The strong elec- 
tric fields induced by the DNA molecule (:^ l.OV/ nm) 
result in Vr 0.2 to 2 meV (or 10~^t), in agreement 
with estimates in the literature [l, [lol, 13 • these 
small Vr values and moderate length, one can then write 
P PmL'^/^0^^ which gives a modest polarization of 
P 0.003% for L = 10. Additional SOI induced by the 
CNT curvature is expected to enhance the spin polariza- 
tion, and may well result in larger P values than these 
minimal estimates suggest. Measuring this effect would 
likely require low temperatures in experiments, although 
increasing the number of wrappings can result in a sub- 
stantial increase in polarization. 
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FIG. 5. (Color online) Normalized polarization P of hybrid 
with (9,0) CNT as function of number of DNA wrappings and 
different Vr for fixed Fermi energy of E=0.03t and Cia = O.lt. 
Solid and dashed lines represent results for wrapping angles 
= 7t/6 and — — tt/S, respectively. Dot-dashed black line 
shows null polarization for the molecule laying parallel to the 
tube. Reversing direction of wrapping reverses conductance 
polarization. Inset shows oscillatory behavior of P for = 
— 7r/3, Vr = 0.2t, for larger range of DNA wrapping. 

In summary, we find spin polarized conductance for 
DNA wrapped CNTs due to the presence of helicoidal 



fields and Rashba spin-orbit coupling-despite the absence 
of magnetic fields. This polarization is sensitive to the 
wrapping pitch and direction, reversing sign for electrons 
and holes; it can be enhanced by orders of magnitude for 
a few wrappings. This remarkable behavior can be used 
in spintronic devices, such as a spin selective filter or in 
identifying the conformation of molecules adsorbed on 
these devices. These results also suggest that other sys- 
tems with hehcoidal Rashba fields would result in conduc- 
tance polarization without applied magnetic fields, a gen- 
eral result that could in principle be implemented in other 
electronic systems, such as semiconductor nanowires. 
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